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ABSTRACT 

Recent large-scale galaxy spectroscopic surveys, such as the Sloan Digital Sky Survey (SDSS), 
enable us to execute a systematic, relatively-unbiased search for galaxy clusters. Such surveys make it 
possible to measure the 3-d distribution of galaxies but are hampered by the incompleteness problem 
due to fiber collisions. In this study we aim to develop a density measuring technique that alleviates 
the problem and derives densities more accurately by adding additional cluster member galaxies that 
follow optical color-magnitude relations for the given redshift. The new density measured with both 
spectroscopic and photometric data shows a good agreement with apparent information on cluster 
images and is supported by follow-up observations. By adopting this new method, a total of 924 
robust galaxy clusters are found from the SDSS DR5 database in the redshift range 0.05 < z < 0.1, 
of which 212 are new. Local maximum-density galaxies successfully represent cluster centers. We 
provide the cluster catalogue including a number of cluster parameters. 
Subject headings: catalogs — surveys — galaxies: clusters: general — galaxies: general 
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1. INTRODUCTION 

The role of environment in the formation and 
evolution of galaxies is a topic o f muc h interest in 
modern astrophysics. iDressled ()1980f ) found that 
the abundance of early-type galaxies increases in 
dense environment indicating that environment has 
an effect on galaxy morphology. This morphology- 
density relat ion has been further studied by many 
authors fe.g.. IWhitmore fc Gilmorelll991l: iDressler et al 



199^ iCoto et al.l 1200.1 iCapak et all 120071 : iPark et al 
2007t l and has been extended to lower density r egions 
and high-rcd shift clusters (iPostman fc Gelled 1198 



iGiovancUi. Havnes. fc Chincarinil 1X9861 : iTullT 

There have been several physical 
posed to explain this relatio ii, 



mechanisms pro- 
such as ram pressure 
stripping (iGunn fc Gotll 119721: iFarouki fc Sha piro 1980; 
Puiita fc Nagashima "'1999'; 'Abadi. Moore, fc Bower 
1999.: Quihs. Moore, fc Bower., 2000,: Chung et al. ,2007,) . 



^^MIE^S'^^l 

:alaxy har assment ("Moore et al.l ll996L Il999f ). starvation 
Larson. T inslcv. fc Caldwell 1980) , and interaction and 
merging (Toomrc fc Toomrc 1972^ Makino fc Hut 1997 



tidal forces (iBvrd fc ValtoneiJ 119901: IValluril 



[Bekki 1999 : Bos elli fc Gavazzi 2006.: Chung et al. 2007 

Some of these processes involve violent gas dy- 
namics, and hence it is quite likely that environ- 
ment influences the star formation ra t e in galax- 
ies (Gislcr 1978; Larson fc Tmslev" 197^; Lewis et alj 
2003; Ch ristlein fc Zabludofi 200i; 

The cosmic star formation his- 
suspected to depend on environ- 
Perhaps 
as 



20021: .Gome z et al. 
Elbaz et all l2007t ). 



tory i tself is also 

ment (jThomas. Maraston. fc Bender! 120051 ). 

as hints of these suspicions, galaxy properties such 



color, lumi nosity, and s i ze app ea r to depend on en- 
vironment (jBalogh et all l2004al |bl: iBlanton et all 120051 : 
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Croton et aDl2005t IPostman et all [20051 : iQuintero et all 
20061 : [Park et al.ll2007| ). However, the exact physical pro- 



cesses that link between galaxy properties and environ- 
ment are still far from being clear. 

Galaxy clusters are the most dense regions in the large- 
scale universe. The number density of galaxies in clusters 
can be several hundred times larger than that of the field. 
They are thus ideal laboratories for studying the effect of 
environment on the formation and evolution of galaxies. 

The first galaxy cluster catalogues (lAbell ,19581 : 
iZwickv et anil961-1968HAbell. Corwin. fc 01owinlll989f l 
were based on the visual inspection of photographic 
plates. While these catalogues have been widely 
adopted, they suffer from many selection effects and 
spurious detections due to projection effects. Fur- 
thermore, the visual inspection of large parts of the 
sky is very time consuming. Systematic searchers 
for galaxy clusters are needed to ensure the efficiency 
and reliability of the cluster detection. Among the 
most popular in the absence of spectroscopic informa- 
tion were the technique using the early- t ype g alaxy 
color-magnitude relation (jGladders fc Ye e! '2000 D and 
findin g brightest cluster ga laxies (BCGs; Bahcal l et al.l 
120031: Ivon der Linden et al.l l2007f). The first automatic 
search f or optical c luster s looking for overdense regions 
was by IShectmanI (Il985fl followed by riumerous others 
Lumsden et al.lll99 2':'Dalto n et allll994j: IPostman et al " 



1996; Croft et al. 



20031 : iMiUer et al 



1997; Kim et al.ll2002l: iKochanek et al 

l200_5) 



With the advent of large-area, spectroscopic redshift 
surveys, such as the Sloan Digital Sk y Survey (SDSS; 
lYork et al.ir2000l : IStoughton et al.ll2002f ) an d Two Degree 
Field Galaxy redshift Survey (2dFGRS; iColless et all 
1200 If ) , large cosmological volumes with redshift have be- 
come available. We are finally able to systematically 
search for galaxy clusters in a large volume by measuring 
the local 3-d number density of galaxies. While the 3-d 
densities measured based on the spectroscopic database 
are more powerful than projection-based 2-d searches for 
delineating the local galaxy distribution, they require a 
high level of coverage of member galaxies in the spectro- 
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scopic survey a priori. However, even the most up-to- 
date surveys (such as the SDSS) show only ~ 60-70% 
completeness rate in dense regions, hampering us from 
measuring the densities accurately. In this paper, we 
introduce an improved method for finding galaxy clus- 
ters by measuring local densities using both spectro- 
scopic and photometric data, and provide a new cata- 
logue of the clusters we have found. We assume cosmo- 
logical parameters flm = 0.3, i^A = 0.7, qo = —0.55, and 
Hq = 70 km s~^ Mpc~^ throughout this paper, and all 
distances are comoving. 

2. THE DATA 

The SDSS is performing a survey to cover a quarter 
of the whole sky. The imaging survey of the SDSS DR5 
contains 215 million objects in 8000 deg^. The SDSS 
spectroscopic survey mapped 5740 deg^ obtaining about 
a million spectra of which 674,749 objects are classified as 
galaxies. The galaxy redshifts provided make it possible 
to study 3-d structure of galaxy distribution. 

The SDSS photometric pipeline provides several kinds 
of magnitudes. For galaxy colors, we use the modelMags 
as they provide unbiased colors regardless of any color 
gradients. For galaxy luminosities, we use petroMags 
as these are a better estimate of the total luminosity 
(jStoughton et al.ll2002D. We also a pply a K-correction as 
described in^ Blanton et aLl([ 2003b') and correct for Galac- 
tic ex tinction with the [Schlcgcl, Finkbcincr, & Davis 
()1998[) values provided by the SDSS pipeline. 

We use the SDSS DR5 both photometric and spec- 
troscopic data. For the galaxies with spectra, we ex- 
tract all types of galaxies in the range 0.05 < z < 0.1 
and having Tpetro < 17.77 (jStrauss et all l2002f ). At 
z = 0.1, r = 11 .n corresponds to an absolute magnitude 
of Mr = —20.55. In order to create a volume-limited 
sample and so avoid biases with redshift, we cut at this 
absolute magnitude. 

For galaxies without spectra, we select all galaxies with 
13.00 < Tpetro < 17.77 to have the same apparent magni- 
tude cut as the spectroscopic data. Its faint limit comes 
from the SDSS spectroscopic survey limit, and the bright 
limit is from the fact that the objects of Tpetro < 13.00 
and z > 0.05 are almost always stars rather than galax- 
ies. 

3. MEASURING GALAXY DENSITY 

3.1. Local Density Measurement 

Popular methods to measure galaxy local densities 
include a simple estimation of spatial number density 
in a certain radius and the distance to the nth near- 
est galaxy. The number density, however, cannot tell 
us about the concentration status which is critical for 
studying dynamical evolution, while using the nth near- 
est galaxy can be biased by local density fluctuations. 
Thus, in addition to the number density, the spatial 
separation between galaxies can be taken into account 
as well to provide additional information. For exam- 
ple, close r neighbors can b e weig hted more than distant 
galaxies. ISchawinski et al.] (f2006) (hereafter S06) intro- 
duced such a weighting scheme with a Gaussian filter. 
The S06 scheme considers neighboring galaxies in the el- 
lipsoid defined by: 



3ct 



3c. (7 



< 1, 



(1) 



where a is the searching distance criterion, is the pro- 
jected distance and rz is the line-of-sight distance, all in 
Mpc, to a neighboring galaxy. It allows all the galax- 
ies within the ellipsoid to be counted for measuring the 
density. The Cz factor is a simple compensation for "the 
finger-of-god" effect due to the peculiar motion and is 
estimated by counting the number of galaxies n within 
the sphere of radius a as: 



1 -|-0.2n 



(2) 



where n is capped at 10. Hence, this can scale along the 
radial direction up to a factor of 3. Density measures 
by adopting a fixed volume in comoving space can lead 
to sampling vastly different volumes in the field and in 
clusters due to the stretching of the line of sight direction 
by the finger-of-god effect. 

The density parameter to each member galaxy based 
on the spectroscopic data is calculated as 
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by summing up for all the members within the search 
ellipsoid. The target galaxy is not included in the density 
summation. Our algorithm estimates the local density 
with surrounding galaxies in a relatively large area and 
thus minimizes local effects. 

In this paper, we attempt to improve the density mea- 
surement scheme of S06 by adding photometric member 
candidates and build on the S06 algorithm to design an 
effective cluster search method. 

3.2. Member Selection with Spectroscopic Data 

Since we do not have the line-of-sight information on 
the photometric member galaxies, we first measure the 2- 
d densities for our spectroscopic galaxy samples ignoring 
the line-of-sight information. 

We first define an ellipsoid as in S06: 



(4) 



We choose ci = 1 Mpc because it optimizes the search 
for the density peak and the brightest cluster galaxy of 
the cluster. For (T2 we adopt the distance corresponding 
to three times the velocity dispersion of the cluster. Note 
that tT2 is used only for selecting member galaxies but not 
for the 2-D density measurement. Readers are referred 
to §3.5 for more details. Galaxies within the ellipsoid are 
regarded as member galaxies of the cluster. Then the 
new 2-d density pspcc for the spectroscopic members is: 
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27rCTi 
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(5) 



where ci = 1 Mpc. 



3.3. Member Selection with Photometric Data 

The SDSS spectroscopic survey tries to cover galax- 
ies as com pletely as possible with a tiling algorithm 
(jBlanton et al. 2003a). However, its selection scheme can 
leave some galaxies unobserved due to fiber collisions. 
Figure [1] shows the extreme case of Abell 2670 which 
shows only 65% of spectroscopic coverage rate within the 
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Fig. 1. — The spectroscopic completeness for the case of Abell 
2670 with rpctro < 17.77. The dots with rectangles are galaxies 
covered by the SDSS spectroscopic survey and the crosses are the 
missed ones. The completeness rate for this cluster is only 55% 
in this field of view and 65% within 1 Mpc radius from the cluster 
center. This figure illustrates the problem of fiber collisions in 
dense environments. 



Fig. 2. — Empirical CMR for the redshift range 0.07 < z < 0.08. 
The gray-solid line shows a linear fit after 2cr clipping until the 
number of galaxies on the CMR becomes constant and the dashed 
lines show the scatter (the gray-solid line ±3(t in g — r color). The 
vertical dot-dashed line is a absolute magnitude cut in our data 
selection. 



1 Mpc radius from the cluster center. This effect is gen- 
erally worse for denser regions because the fiber collision 
problem is obviously worse in the more crowded regions. 
We estimate the spectroscopic completeness, that is, the 
number fraction of galaxies covered by the SDSS DR5 
spectroscopic survey, to be /spec ~ 65% for rich clusters. 
An independent study has also reported that 30% of true 
brightes t cluster galaxies are missed by the spectroscopic 
survey (|von der Linden et al.l 120071 ). The "incomplete- 
ness" of the spectroscopic survey causes a problem in 
measuring the densities of galaxies in dense regions. We 
attempt to alleviate this by further considering photo- 
metric data. 

Early-type galaxies in clusters have a tight correla- 
tion between their optical colors and luminosity, known 
as the color- mag nitude r elatio n (CMR). The CMR was 
first observed bv iBaunj (I1959D and is often reported to 
be universa.1 with i n errors (IVisvanathan fc Sand agdl 19771: 
iHogget al.l[2004t iLopez-Cruz. Barkhouse-TY ec 200"^. 
Based on these assumptions, we can select candidate 
cluster members even though they do not have redshift 
because galaxies on the CMR have a high likelihood of 
belonging to the cluster. Thus, we find the slope and 
scatter of the CMR in order to select "photometric mem- 
bers". For small clusters, the slope of the CMR can be 
difficult to achieve. Hence, we construct an empirical 
CMR by stacking galaxies that are members of 20 typ- 
ical clusters determined by their redshifts. We do this 
for five redshift bins with Az = 0.01, and an example is 
shown in Figure [2] A linear fit with 2 a clipping is iter- 
ated until the number of galaxies on the CMR remains 
constant. The galaxies residing on the CMR ±(3 x rms) 
in g — r color are selected as "photometric members" . 
We compute the 2-d photometric density Pphot following 
Eq. [5] using only the photometric members. 

The total density p of a galaxy is calculated by com- 
bining the spectroscopic and photometric density param- 
eters: 

P = Pspcc + /Ophot- (6) 



In the process of combining pspcc and pphot we are los- 
ing the line-of-sight information and our density mea- 
sures are projected into only two dimensions; but, our 
method still provides improved density measures as will 
be demonstrated in §4. through our follow-up spectro- 
scopic observations. 

When we determine whether a particular galaxy with- 
out redshift belongs to a cluster as described in this sec- 
tion, we assume that the galaxy is at the redshift of the 
cluster and apply the same absolute magnitude cut as 
discussed in §2. 

We have attempted improving our local density mea- 
sure further by applying a galaxy luminosity weight (i.e. 
weighting more luminous galaxies more) but found no 
significant difference. Hence, we have decided to ignore 
it. 

3.4. The Efficiency of the CMR technique 

We expect that the CMR to be effective for finding 
cluster member galaxies but that it may still suffer from 
the projection effect of background and foreground ob- 
jects. Hence, we test its efficiency. We restrict the test 
on the spectroscopic data only, because only with red- 
shift we can determine whether a particular galaxy is a 
cluster member or not, that is, for a convincingly high 
confidence. Besides, we use early-type galaxies alone for 
this test as our CMR method works only for them. In or- 
der to select early- type galaxies (ETGs) , we use the mor- 
phological index fracDev from the SDSS pipeline, which 
shows the weight of the de Vaucouleurs profile in the two- 
component profile fit with the de Vaucouleurs and expo- 
nential profiles. We apply the highly conservative r-band 
fracDev > 0.95 cut, and the galaxies selected are highly 
likely to be early type. Using the spectroscopic members 
meeting the same fracDev criterion, we check the valid- 
ity of our scheme defining the following two CMR-related 
quantities. 



Completeness - 



# of ETC members in CMR 
# of all ETC cluster members 



X IQCF) 
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^ . # of ETG members in CMR 
P""^y= # of ETG galaxies in CMR ^ ^^^^^^ 

The "CMR completeness" is a measure of the num- 
ber of the early-type cluster members that also sit on 
the CMR. The "CMR purity" tells us the fraction of 
galaxies selected by the CMR method that are indeed 
early-type cluster members. Figure [3] presents the CMR 
completeness and purity of our clusters as a function of 
density. The CMR completeness is around 90% mean- 
ing that 90% of the early-type (by fracDev) member 
galaxies indeed sit on the CMR. The purity is 70-85% 
meaning that only 70-85% of the galaxies on the CMR 
are early-type (by fracDev) members. The rest may 
be background or foreground galaxies. Exact values of 
these parameters modestly depend on the fracDev cri- 
terion adopted, but on the whole both parameters show 
high values validating our scheme. 

3.5. Finding Clusters 

The SDSS DR5 photometric data includes about a mil- 
lion galaxies in the magnitude range described in §2. 
Searching for clusters through the entire data set would 
require an exhaustive amount of effort. So we begin 
our cluster search first by estimating the local density of 
spectroscopic members and identify about 7000 galax- 
ies with high density measures as candidate cluster posi- 
tions. For this selection we use 3-D spectroscopic density 
parameter (pspec.sD > 0.5). 

We then select all the photometric and spectroscopic 
galaxies within a 2 Mpc radius from the cluster center 
candidate galaxies and estimate their total 2-d local den- 
sity with both the spectroscopic and photometric mem- 
bers for all galaxies in this radius. 

With these 2-d spectro-photometric density measure- 
ments, we find the galaxy with the highest density in 
the area of a 2 Mpc circle and Az = 0.01. The highest 
density galaxy in a local area is defined as the maximum- 
density galaxy (hereafter MDG). As the MDG resides at 
the densest part of the selected volume, we define it as 
the cluster center and declare it a cluster if the p of an 
MDG is greater than 4.0 which roughly corresponds to a 
cluster velocity dispersion of 200kms~^. 

We expect that the MDG should be close to the 
cluster center an d coincide with the b rightest clus- 
ter galaxy (BCG; [Matthews. Morgan, fc Schmidt Il964l : 
lOegerle fc Hillll2001h . We compared the projected sep- 
aration of BCGs to MDGs and found that they do not 
always coincide, however. We assume that a good clus- 
ter search scheme should identify MDGs as BCGs as well, 
and thus we vary the member search radius (the short 
axis of the ellipsoid) hoping to minimize the separation. 

Figure 2] shows the separation between MDGs and 
BCGs. The mean and error bars reflect the large scat- 
ter in the sample. The median is a better estimate of 
the typical distance and shows a convergence towards a 
lower value of search radius. But too small a value of the 
search radius would miss too many member galaxies in 
the density measurement and be dominated by local den- 
sity fluctuations. The typical Virial radius of our clusters 
is 0.5 Mpc in the 2-d projection, and a substantial frac- 
tion {r-^ 50%) of clusters have a Virial radius greater than 
0.5 Mpc and extend beyond IMpc. Considering this, we 
have chosen 1 Mpc as our search radius. It should be 
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Fig. 3. — The solid and dashed hnes show completeness and 
purity of the photometric member candidates found from our tech- 
nique (see the text for definition). The errors are from the Pois- 
son statistics in the equal number bin. The rectangles are shifted 
slightly to the left for the clarity of the figure. The completeness 
is constantly ~ 90% and the purity is 70-85% with an apparent 
density dependence. The high values of these parameters validate 
our CMR technique. 
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Fig. 4. — The projected separation between brightest cluster 
galaxies (BCGs) and maximum-density galaxies (MDGs). Trian- 
gles are the mean and squares are the median separations. The 
error bars show the standard deviation for the mean and the me- 
dian absolute deviation for the median. When the search radius to 
find member galaxies is 1 Mpc, the separation is small enough that 
the cluster search is effective. But a substantially smaller values 
than 1 Mpc would miss many memeber galaxies because the typical 
size of clusters can be easily as large as 1 Mpc (see the text). 



noted that our search using spectroscopic data and the 
CMR is liable to missing clusters that are dominated by 
blue, starforming galaxies, if any. 

4. RESULTS 
4.1. Cluster Catalogue and Properties 

We provide our cluster catalogue in Table [H We have 
found 924 galaxy cluster candidates with p > 4.0 which 
roughly corresponds to 200 km s^^. Among them, 212 
clusters are newly found in this study. 

Figure [5] shows a sample from the newly-found clusters. 
The BCG and MDG are marked as number 1 and 7, 
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TABLE 1 
Cluster catalogue with p > 4 
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1 


258.20851 


64.05295 


0.08257 


32.197 


148 


2.785 


1299 


258.12006 


64.06076 


0.07340 


A2255 


2 


358.55746 


-10.39042 


0.07623 


29.234 


91 


1.842 


840 


358.55703 


-10.41904 


0.07766 


A2670 


3 
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5.14974 
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A1809 


4 
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1.487 
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4.38621 


0.00000 


A2048 
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255.66908 


33.52032 
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22.783 
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1.813 


1135 


255.63809 


33.51666 


0.08640 


A2245 
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27.23342 


0.08969 


22.476 
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239.58334 
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A2142 
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0.07688 
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205.54018 
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0.00000 
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186.91015 
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41 
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230.33576 


30.67093 
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0.07845 
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10 


255.64172 


34.07809 


0.09871 


19.178 


64 


1.546 


1120 


255.67708 


34.06003 


0.09891 


A2244 



Note. — Column 1: ID; Column 2,3: the coordinates of the maximum-density galaxies (MDGs); Column 4: the redshift of clusters 
estimated with spectroscopic data; Column 5: the density of the cluster MDG within 1 Mpc; Column 6: the number of member galaxies 
of clusters within -R200; Column 7: the virial radius of clusters obtained by the number density profile of the galaxy clusters; Column 8: 
the velocity dispersion within -R200; Column 9,10: the coordinates of the brightest-cluster galaxies (BCGs); Column 11: the redshift of the 
BCGs. The entries with zero mean the absence of the redshift information; Column 12: matching result with known catalogue of galaxy 
clusters. The complete version of the galaxy clusters is available on the web-page |http:/ /gem.yonsei.ac.kr/html/cluster.php| 
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Fig. 5. — A sample rich cluster that is newly found in our study. The density ranking is 142th, out of 924, with the density parameter 
p = 9.1 and Afgai = 16 within 1 Mpc, where its virial radius -R200 is estimated to be 0.737Mpc and N200 = H. The size of each circle 
denotes its r-band luminosity. Blue and red circles show photometric and spectroscopic members, respectively. SDSS image cut-outs are 
shown for the 16 members and the central MDG (number 7). The BCG is marked as number 1. 
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respectively. Blue and red circles denote photometric 
and spectroscopic members, while the size of circles is 
proportional to the galaxy r-band brightness. Postage 
images of the member galaxies are also shown. Its MDG 
is ranked 142th out of 924, with p = 9.1 and iVgai = 16 
within 1 Mpc. For comparison, the central galaxy in the 
Virgo cluster, M87, has p = 7.5 and A'gai = 19 within 
1 Mpc, and thus is comparable. 

Table [1] provides various properties of our galaxy clus- 
ters. The positions of MDGs and BCGs are included. 
The redshift of each cluster is obtained by finding the 
redshift peak of the spectroscopic members and fitting 
it by Gaussian function. When the number of spectro- 
scopic members is less than 4, we simply find median 
redshift of spectroscopic member galaxies. We calcu- 
late the virial radius r2oo by counting the number den- 
sity of member galaxies and comparing it to the cos- 
mological background number density of galaxies. We 
estimate the number of galaxies within a virial radius 
-^200 by counting the galaxies within r2oo- The ve- 
locity dispersion is calcula.ted w ith biweight estimator 
(jBeers. Flvnn. fc Gebhard^ 119901 ) in r2oo after 3-sigma 
clipping. As a proxy to the cluster size and density, we 
can use p regardless of the cluster size (virial radius). 

In §3.5, we introduced the separation between the 
MDGs and BCGs as a test of how well the cluster search 
scheme finds cluster centers. In Figure [6l we show the 
separation between MDGs and BCGs in our method with 
(bottom) and without (top) added photometric mem- 
bers. It illustrates the effectiveness of the added photo- 
metric members. 

The galaxy number density profile is a good 
proxy for the mass profile a nd thus dark mat- 
ter profile of galaxy cl u sters (ICarlberg et al.l IT997I: 
van der Marel et al. '2000'; 'Katgert, Biviano, & Mazur^ 
2004; Lin, Mohr, & Stanford 2004; Hansen ct al. 200^ 
In Figure [71 we show the galaxy number density 
profile of two examples. We compare them with 
the Navarro- Frenk- White (NFW ) dark matter profile 

avarro. Frenk. fc Whitelll996l ). which is defined as 

= (9) 
Pc {r/rs){l + r/rsY'' 

and can be projected into a 2-dimensional space as 



= 2710 



5r. 



;de (10) 



Pc Jq {rp/rsCos9){l + rp/ cos6 

where uq is normalization factor, rp is projected distance 
and 9 = cos^^-^. The number density profile is in rea- 

'"200 ^ 

sonablc agreement with the NFW profile as found also by 
[Hansen et al. (2005). The galaxy number density profile 
only based on the spectroscopic members is lower than 
the one with photometric members added. The depar- 
ture gets larger in the central region of the clusters where 
the spectroscopic coverage is poor. 

A larger cluster would have a higher value of p, the 
density measure within 1 Mpc. We show the distribution 
of p, N20Q and the velocity dispersion of clusters within 
r2oo in Figure[8l The (minimum, median, and maximum) 
for the four panels are (0, 5, 190) for A^200, (4.0 , 5.7, 32.2) 
for p, (0, 263, 1299) for u, and the redshift range is 0.05 
through 0.1. The density and the velo city dispersion 
of M87 (.Binggeh. Tammann fc Sandagel[l987.) in Virgo 
cluster is shown as well. 
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Fig. 6. — We show the separation between BCGs and MDGs. 
The dotted Unes represent the median separation. The top panel 
shows the separation when the CMR-selected photometric mem- 
bers are not included in the density measurement. When we 
add the photometric members (bottom), the separation gets much 
smaller yielding a more reliable cluster center detection. 
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Fig. 7. — The projected radial profiles of the galaxy number 
density of sample rich clusters centered at their BCGs. The plus 
signs show the galaxy number density profiles from our scheme 
(including both spectroscopic and photometric members), and the 
triangles are those with only spectroscopic members. Also shown 
for comparison are the projected NFW halo profiles fitted to our 
measurements with both spectroscopic and photometric members. 
Our measurements are in good agreement with the projected NFW 
profiles. 
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Fig. 8. — Histogram of N200, P, velocity dispersion, and redshift. We define galaxy clusters when a condition p > 4.0 is met. In the top 
left panel two clusters of N200 = 190 and 148 are out of bound. The vertical dotted lines show median values. The density and velocity 
dispersion of M87 are shown as the arrow. 



In Figure [5] we compare p with iV200- Also shown is 
the linear fit to the data. The density parameter p ap- 
pears to scale as a power law with N200 rather than with 
a linear relation. That is, in dense regions, the linear fit 
would have a steeper slope. This is because the Gaussian 
weight in our density measuring scheme takes the dis- 
tance to each member into account. For a given number 
of cluster members, a more concentrated cluster would 
have a higher value of p. 

In Figure [To] we compare between various cluster prop- 
erties. The density parameter is a good proxy to the size 
of galaxy clusters. It also shows a good correlation with 
the velocity dispersion of clusters. Thus, our density pa- 
rameter can be an obvious indicator of dynamical mass 
of a galaxy cluster as well, as shown in the bottom panel. 

4.2. Comparison to other catalogue 

We perform a cross-matching between existing cata- 
logues and ours. The matching criteria are 2 Mpc in the 
projected radius and Az = 0.01. When the known clus- 
ters do not have redshift (e.g., Zwicky clusters), we just 
match them within 2 Mpc radius in the sky. Each cluster 



TABLE 2 

The number of galaxy clusters found in our scheme 
matched to previous catalogues. 



Cluster Catalogue 


Number 


Recovery (%) 


AbcU 


128 


79% 


Zwicky 


457 




C4 


300 


66% 


Others 


889 




Newly Found 


212 




Total 


924 





Note. — The existing cluster catalogues are from the VizieR 
Catalogue Service. 

in our catalogue is matched with other catalogues sepa- 
rately. The clusters that are detected in our method are 
generally in the existing catalogues (See Table [2]). Most 
rich clusters are matched with Abell or Zwicky clusters 
and this implies that the new algorithm using MDGs 
successfully finds rich galaxy clusters. For those clus- 
ters in the existing catalogues which include redshift, we 
can check the rate at which our method recovers clusters 
known previously. Some Abell clusters have redshifts and 
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Fig. 9. — The density parameter p is proportional to N200 
scaling as power law. Note that the density parameter includes 
the distance information in addition to the number. The thick line 
is the linear fit to the data and the grey band shows the standard 
error of the fit. 




0.6 0.8 1.0 1.2 1.4 

log p 



Fig. 10. — The density parameter p is in agreement with the 
properties representing the mass and size of clusters. The virial 
radius in the top panels has a tight correlation with p. The density 
parameter correlates with the velocity dispersion and dynamical 
mass in the middle and bottom panels. The thick gray line is our 
linear fit and dy/dx presents the exponents of the power law. The 
error bars represent mean scatter. 



the C4 catalogue clusters always do (jMiller et al.ll2005[) . 
We cross-check our clusters against the Abell and C4 
clusters that are in the SDSS area and in our redshift 
range. 128 out of 162 Abell clusters with redshifts and 
300 out of 458 C4 clusters are detected in our method, 
corresponding to detection rates of 79% and 66%, respec- 
tively. If we lower the p cut for "clusters" from 4.0 to 
1.0 we would find 89% and 95% of Abell and C4 clus- 
ters, respectively. For example, by using a, p — 1 cut, we 
can recover additional 18 Abell clusters. However, such 
a low value of cut {p ^ 1.0) would also find significantly 
more cluster candidates (3,959) many of which may be 
false detection. Our density measures with p = 4 cut are 
in general based on about 7 member galaxies. Besides, 
when only a small number of galaxies are used for density 
measurement, the density measures become less reliable, 
too. Since our goal is to construct a database of galaxy 
clusters with improved and hence reliable density mea- 
sures, we think our p = 4.0 cut is a good compromise. 
However, we admit that this choice is still somewhat ar- 
bitrary and by using this cut we are neglecting smaller 
clusters. In this regard, our 924 cluster candidates rep- 
resent a robust sample of relatively-rich clusters. 

The galaxy clusters that were known to Abell but 
missed by our method are mostly due to the p cut as 
mentioned above. In addition, our method missed 9 
Abell clusters which are mainly composed of relatively 
faint galaxies that do not make it to our volume-limited 
sample. Two Abell clusters were missed because they are 
found by our method to be parts of larger clusters, just 
like the M87 and M49 groups in Virgo cluster. Five addi- 
tional clusters were missed due to various minor reasons: 
e.g., mismatch in the redshift of the center galaxy of a 
cluster, mismatch in the position of the center galaxy of 
a cluster, and so on. 

We c ompare our new ca talogue to the C4 cluster cat- 
alogue (IMiller et all 120051) . While we use the density p 
to represent the richness of a cluster, the C4 Catalogue 
provides the number of galaxies within various radii. We 
compare our density measures within 1 Mpc, p, with the 
C4 richness (C4 parameter wmaglOOO) in Figure [TTl The 
match is reasonably good but shows a large scatter. The 
systematic difference in vertical scale has a couple of crit- 
ical origins. Most importantly, our limiting magnitude, 
Afj. ^ —20.5, is somewhat brighter than that of C4 (- 
19.9). Secondly, our density is not a simple count of 
galaxies but weighted by the distance from the cluster 
center. Thirdly, our method, unlike C4, includes addi- 
tional photometric member candidates. Fourthly, our 
choice for the search radius (1 Mpc) is smaller than that 
of C4 (l//iMpc). Lastly, the member-search volume in 
C4 is a cylinder with a line-of-sight length correspond- 
ing to four times the velocity dispersion, while ours is 
an ellipsoid with a length that is 3 times the velocity 
dispersion. The large scatter is mainly a result of two 
effects. To begin with, the position of the central galaxy 
of a cluster can be different mainly because our method 
includes photometric member candidates. In addition, 
the redshift of a cluster can also be different. We notice 
some difference in the velocity dispersion measured as 
well, which also contributes to the scatter. 



4.3. Effects of photometric members 
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Fig. 11. — Comparison between catalogues in terms of the 
richness of the clusters. X-axis shows our density parameter p 
(listed in Table 1) and y-axis shows Ngal (wmag 1000) from the 
C4 catalogue. The large scatter has several origins (see text for 
details). The line shows the linear fit. 



The inclusion of the CMR-selected photometric mem- 
bers increases our local density measures by ~ 22% (See 
Figure [T^ . In this figure, the small black dots show the 
density measures of our clusters based only on the SDSS 
database. A linear fit to the data is shown as a dashed 
line (showing the 22% increase) and the one-to-one rela- 
tion is presented as a dotted line. 

To check the validity of the CMR technique further, 
we attempt to recover our spectro-photometric den- 
sity approximately by dividing the measured spectro- 
scopic density pspcc by the estimated spectroscopic com- 
pleteness fraction /spec- Figure [T3l compares Pspoc//spcc 
(crosses) with our spectro-photometric density p (the up- 
per ends of arrows). They agree reasonably well. Since p 
is a complex parameter including the distance informa- 
tion while /spec is not, we also make the comparison to 
A'gai, which also shows a tight correspondence. 

4.4. Additional spectroscopy for the test 

For 22 of our clusters we were able to find 2dF spec- 
troscopic data with which we identified further member 
galaxies. Using these new spectroscopic members, we 
computed new (enhanced) Pspcc and compared them to 
our SDSS-spectro-photometric p in Figure fTSl The x axis 
shows our 2-d SDSS-spectroscopic density (pspcc) that 
suffers from the incompleteness problem. The y axis 
shows density measures improved upon pspoc through 
various methods. The blue filled circles show the new 
spectroscopic density measures for the 22 clusters includ- 
ing the added 2dF member candidates: Pspcc.SDSS+2dF- 
Only a few of the new SDSS-|-2dF spectroscopic den- 
sity measures (pspoc,SDSS+2dF) overshoot our 22% en- 
hancement prediction (dashed line). Even in this case, 
however, the SDSS-|-2dF spectroscopic density mea- 
sures (pspec,SDSS-i-2dF) are Still lower than our SDSS-only 
spectro-photometric density measures (pspcc+phot.SDSs; 
that is, our density measures in Eq. 6) in most cases, as 
indicated by the arrows. The end points of the arrows in- 
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Fig. 12. — We compare our density measure p with pspcc (the 
density only with spectroscopic data). The dotted line presents the 
one-to-one relation between p and pspoc and the dashed line shows 
a linear fit with the y-intercept fixed at the origin. The departure 
gets larger in dense regions. The density increases by 22% through 
the addition of photometric member galaxies. 



dicate our SDSS-spectro-photometric density measures. 

In order to test our result, we have performed CTIO 
Hydra MOS observation of Abell 2670 and significantly 
increased the spectroscopic coverage of this cluster from 
65% (within IMpc) to 92%. We chose this cluster as 
a test case because it suffers from an extreme case of 
the spectroscopic incompleteness problem as mentioned 
in §1. 

The observation was performed on 1-3 December 
2006 with Hydra, a Multi-Object-Spectrograph mounted 
on Blanco 4m- Telescope at CTIO. The wavelength 
range covers 3600-8000A and the spectral resolution is 
2.3A/pixel with mean signal-to noise ratio of 13. Orig- 
inally all of the early-type galaxies in Abell 2670 with 
r < 18 were selected as targets. In order to execute 
our plan as completely as possible, we used three fiber- 
configurations to evade fiber collision issues in dense re- 
gions. For each fiber configuration, a total of 45 min- 
utes which consists of 3 exposures of 15 minutes was 
taken for median combine and cosmic-ray removal. All 
the processes are done mainly with IRAF hydra pack- 
age. The radial velocities of galaxies were measured by 
cross-correlation method using IRAF rvsao package and 
cross-correlation templates of SDSS. In the top right of 
Fig- [El we show our new SDSS-I-CTIO spectroscopic 
density (pspec,SDSS-i-CTio) at w 25 (purple diamond at 
the top right). The value, based on a relatively com- 
plete (92%) spectroscopic survey, now is a factor of two 
larger than the previous Pspec.SDSS and much closer to 
our spectro-photometric density Pspcc+phot.SDSS- 

5. SUMMARY 

We develop and test a method for finding galaxy 
clusters in the SDSS spectroscopic and photometric 
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database. Our method improves over the previous den- 
sity measurements by finding additional cluster member 
candidates via the color-magnitude relation technique. 
The problem of spectroscopic incompleteness due to fiber 
collisions in dense fields is minimized in our method. The 
member galaxies selected by the color-magnitude relation 
lead to a satisfactory completeness and purity. With this 
method, we find 924 galaxy clusters of which 212 are 
new. We provide a catalogue of these clusters including 
important properties such as the virial radius, velocity 
dispersion, and richness parameters. The density we es- 
timate is in good agreement with the properties relating 
to the cluster mass and size. We also provide an estimat- 
ing scheme for our spectro-photometric density measure 
simply using the spectroscopic completeness rate. Our 
new density information on galaxy clusters could be use- 
ful for the study of the environmental effect on the galaxy 
evolution. 
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Fig. 13. — Top: Comparison between the spectroscopy-only 
density measures (x axis) and various improved density measures 
(y axis). The solid triangles show our spectro-photometric densities 
(p). The solid circles show the spectroscopic densities based on the 
SDSS and 2dF database combined. The purple diamond in the 
top right shows the spectroscopic density of Abell2670 based on 
both the SDSS and our new CTIO database. The crosses show the 
result of rough correction for the incompleteness problem, using 
Pspec (i.e., X axis) and the completeness fraction /spec- Bottom: 
the same as the top panel but for the number of galaxies. The 
dotted and dashed lines are the same as in Fig. 12. 
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